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Principal properties of the velocity distribution of dark 
matter particles on the outskirts of the Solar System 
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ABSTRACT 

The velocity distribution of the dark matter particles on the outskirts of the Solar 
System remains unclear. We suggest to determine it using experimentally found prop- 
erties of the oldest halo objects. Indeed, the oldest halo stars and globular clusters 
form a collisionless system, as well as dark matter particles do, and they evolved in 
the same gravitational field. If we accept this analogy, we can show that the velocity 
distribution of the dark matter particles should be highly anisotropic and have a sharp 
maximum near v ^ 500 km/s. The distribution is totally different from the Maxwell 
one. 

We analyze the influence of the distribution function on the results of dark matter 
detection experiments. It is found that the direct detection signal should differ no- 
ticeably from the one calculated from the Maxwell distribution with (v) ~ 220 km/s, 
which is conventional for direct detection experiments (the ratio depends on the de- 
tector properties and typically falls within the range between 6 and 0.2). Moreover, 
the sharp distinction from the Maxwell distribution can be very essential to the ob- 
servations of dark matter annihilation. 

Key words: cosmology: dark matter, elementary particles, methods: analytical. 



1 INTRODUCTION 

One of the most evident manifestations of the dark 
matter existence is the detection of huge invisible ha- 
los (with density profiles p ~ f^^) surrounding galaxies 
iMarochnik fc Suchkov|[T99^ '). Our Galaxy also has such a 
halo. We symbolize the orbital radius of the Solar System, 
the average velocity of the Galaxy rotation, and the escape 
velocity at this radius by tq, vq, and Veac, respectively. We 
also denote the radial and tangential components of a dark 

matter particle velocity by Vr and Vp e 



'-'^ + Vg. Velocity 

distribution of the dark matter particles (hereafter DMPs) 
inside the halo is poorly known; it is us ually supposed to 
be M axwell with a cut-off when v > Veac (|Cerdeno fc Greenl 
I2OIOI ). 
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We accept r© = 8 kpc, vq = 220 km/s, v^sc ~ 643 km/s. 
A'' is a normalizing constant, and for the chosen parameters 
A'' ~ 1.001. It is worthy of noting that Vorb remains almost 
constant throughout the halo. 

Distribution Q faces with difficulties. In fact, in the 
framework of collisionless dynamics it can be naturally ob- 
tained from profile p ~ r~^, subject to the condition, how- 



ever, that function / is isotropic, i.e. / depends only on |ti| 
(so called isothermal model). This assumption seems highly 
improbable. Indeed, if ([l| is true, the majority of the dark 
matter particles has large specific angular momentum p = 
DJl/m^ — [vxr]. The average angular momentum of the par- 
ticles (and, consequently, of the halo) is zero. However, the 
root-mean-square momentum is \/{p?) — 1800 kpc • km/s 
at tq. Moreover, since in this model p ~ rVorb and Vorb is 
constant in the halo if p ~ , the root-mean-square an- 
gular momentum reaches an incredibly huge value \/{p?) ~ 
4 • 10* kpc • km/s at the edge of the halo (r ~ 200 kpc). 
Meanwhile, according to modern cosmological conceptions, 
not only the total angular momentum of the halo but also 
the momentum of each particle should have been negligi- 
bly small on the linear sta ge of the structure formation 
(IGorbunov fc RubakovllioTTI ). The halo could gain some an- 
gular momentum later, as a result of tidal perturbations or 
merging of smaller halos; nu merical simulations show, how- 
ever, that it cannot be large jMaccio et al.|[2007l ). 

Some results of stellar dynamics are frequently used in 
order to show that the dark matter particles could gain a 
large angular momentum during the Galaxy evolution. The 
parallels between DMP and stellar dynamics, however, are 
not universally true. The point is that stars are compact 
objects, and their gravitational field can be strong, at least. 
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locally. On the contrary, the small- scale gravitatio nal field 
of the dark matter is always small (jBaushevI |2009| ) . There- 
fore important relaxation mechanisms of stellar systems like 
close pair approaches or an interaction with the interstellar 
medium are co mpletely ineffective for DMPs. So called vio- 
lent relaxation (jLvnden-Belll [l967l ) is perhaps the only way 
to impart significant angular momentum to the dark matter 
particles. However, it acts also on the halo stars; moreover, 
its efficiency decreases with radius, and it should stronger 
affect the star distribution since the stellar halo is more com- 
pact. 



To summarize: all halo objects initially had Vp 







and later gained some angular momentum because of vari- 
ous processes like relaxation or tidal effects. All the mech- 
anisms increased (Up) of the halo stars, at least, as much 
as of the DMPs, while some mechanisms affected only the 
stars, and not the DMPs. Consequently, velocity distribu- 
tion of the oldest stellar halo population should be closer to 
the Maxwellian one, than the distribution of the dark mat- 
ter particles. In particular, the tangential velocity dispersion 
cr(up) of the DMPs cannot be larger than that of the halo 
stars at the same radius. 

Modern observations of the oldest halo stars - subd - 
warfs - confirm the above reasoning dSmithelaDliooi). 
Their tangential dispersion o{vp) = ctq — 80 km/s, which 
corresponds to ^/JJfi) ~ 900 kpc • km/s, is two times lower 
than in Moreover, the distribution widely differs from 
the Maxwellian: a{vr) is much larger than cr(wp). Conse- 
quently, cr{vp) of the dark matter particles on the outskirts 
of the Solar System does not exceed ao = 80 km/s and can 
be even smaller. Second, the observations show that the halo 
stars have not yet relaxed and their orbits are rather pro- 
late. On the other hand, if distribution ([T]) is correct, the 
ellipticity of the majority of DMP orbits is small, and then 
the dark matter is the only class of halo objects that move 
almost circularly. It seems much more natural to assume the 
opposite, and we come to the premises we use throughout 
this article: 

1) The specific angular momentum ^ of, at least, the 
main part of the particles is fairly small, and their orbits 
are rather prolate (below we will express this supposition 
quantitatively) . 

2) The Galactic halo is stationary and spherically sym- 
metrical. The latter supposition is not quite acc urate: a part 
of da rk matter can form a so-called thick disk (|Read et al.l 
I2OO8I ). moreover, the influence of the star disk also takes 
place. However, our assumption is quite acceptable for an 
estimative consideration. 

3) In accordance with observations (jPutton et al.l 

we assume that the density profile of the Galaxy is 
p oc up to some large enough radius R. The profile can- 
not be valid for an arbitrary large r, since in the opposite 
case the halo mass would be infinite. Starting from some 
radius (we denote it by R), the halo density drops much 
faster. As we will see, this is the main parameter defining 
the DMP velocity distribution. The profile p ~ has an 
another demerit: the predicted annihilation signal (which is 
proportional to p^) diverges when r — > 0. However, as it will 
be shown below, the velocity distribution of DMPs on the 
outskirts of the Solar System weakly depends on the density 
profile inside the solar orbit. Therefore we will not discuss 
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Figure 1. Anisotropic I I13I I (solid line) and Maxwell (([TJ, dashed 
line) distributions of dark matter partieles over the velocity ab- 
solute value. The dotted line represents the velocity distribution 
if the halo edge is smoothed according to II15I I. 



a very complex question of the dark matter profile near the 
galactic centre. 

4) We assume that the dark matter mass outside of 
the radius R is negligibly small as compared with the total 
halo mass. This supposition is the most discussable; how- 
ever, there are stro ng arguments in favour of it . First, exper- 
imental data show l|Binnev fc Tremaindliooj ) that p ~ r 
up to very large distances. Second, in order that the total 
halo mass is finite, the density must fall much faster than 

(at least, faster than r~'^) at large distances. For in- 
stance, one of the most popular Einasto profil e predicts an 
exponential density decreasing (Ein asta Il965l l. Third, the 
outer regions of the halo are respectively weakly bound in 
the gravitational field and can easily be torn away by tidal 
effects. 

These four suppositions turn out to be sufficient to de- 
rive the velocity distribution more or less unambiguously. All 
our calculations are, of course, focused on determination of 
the DMP velocity distribution on the outskirts of the Solar 
System. 



2 CALCULATIONS 

For our Galaxy we accept i? = 210 kpc, which corresponds 
to the total mass of the Galaxy M = 2.3 • 10^^ M©. The 
module of gravitational potential on the edge of the halo 
is equal $ = GM/R. It is easy to see that $ = uq. Since 
p ~ the mass inside some radius can be found as -^Af, 
gravitational field is equal g — G^, and we obtain the 
gravitational potential inside the halo: 
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R 
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Let us start our consideration from the case when the DMPs 
have no angular momentum at all. Then their trajectories 
are radial (w = |wr|), and the task becomes one-dimensional. 
Therefore the particle distribution in the halo can be entirely 
described by a single function ^p{r,v), so that tp{r,v)drdv 
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gives the total mass of dark matter in the element of phase 
space drdv. ip differs from the standard distribution func- 
tion only by a insignificant multiplier — the DMP mass 
rrix- Since we suppose that the dark matter particles are 
propelled only by the gravitational force, rrix is not impor- 
tant for us, the calculations are formally valid even for the 
dark matter consisting of primordial black holes, and such a 
definition of tp allows us to avoid of the undesirable depen- 
dence on the DMP mass. The halo density is bound with 
the function tp hy a trivial relation 



77 = 47rr p{r) — / tp{r,v)dv 
Jo 



(3) 



Here we introduced a more convenient variable rj instead of 
p, ri — const if p ~ r~^. In our case 77 = M/R, and we obtain 
a determining condition on 



il}{r,v)dv = — 
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Each dark matter particle executes a radial oscillation 
around the galactic centre. We denote by ro the maximum 
distance it moves off the centre. Its velocity in potential ((2]) 
is equal to 



\vr\ = 1) = ^/2$ln — 
r 



(5) 



and we obtain useful equations: 
dv 



'-o = rexpl — 



Here T is the time required for the particle to fall from 
ro to the centre. We introduce a distribution function ^ of 
the particles throughout parameter ro, so that ^(ro)dro is 
the total mass of DMPs which apoapsis lies in the interval 
[ro;ro -I- dro]. The r-coordinate of these particles varies be- 
tween and ro, and they give a yield into the halo density 
over all this interval. Indeed, the fraction of time the DMPs 
from the subsystem under consideration pass in an interval 
[r;r -I- dr] is equal to dt/T = dr/{vT). Since the total par- 
ticle mass of the subsystem is 5(ro)dro, the contribution to 
the halo mass in interval dr is equal to 



dM = ^^drdro 
vT 



(7) 



Velocity interval dv that is covered by the particles of the 
subsystem at radius r is rfw = dro. Substituting this to 
(O and taking into account that dM = tp{r,v)drdv, we ob- 
tain the general equation for xp{r,v): 



ip{r,v) = 



vT 



Now we substitute here equations (O and 



i>{r,v) = 



rexp — 



(8) 



(9) 



We can easily find function ^ if we suppose that the halo 
boundary is sharp, i.e. the density obeys the law p ~ 
up to radius R and is equal to zero just after it. Then the 
maximum velocity the halo particles may have at radius r 
is 
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It is a matter of direct verification to prove that function 

M 
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satisfies condition Q. This function has a peculiarity at 
ro = R. So the main part of DMPs comes to us from the 
very edge of the halo. The distribution function is 

2Af 
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where Vmax is defined by (fTO)) . v £ [0;umai]. The dis- 
tribution through the radial velocity Vr can be easily ob- 
tained from (|12p : tp{r,Vr) = ^{r,—Vr) = ^p{r,v)/2, where 
Vr G [—Vmax', Vmax]- Normalized velocity distribution on the 
outskirts of the Solar System is equal to 



~^ 7r^(2.2i;0)2 -t;2 



(13) 



The distribution has a peculiarity at Vmax — 2.55w0 ~ 
562 km/s, which is a result of the supposition that the den- 
sity after r = R immediately drops to zero. In actuality there 
is a characteristic length I of density decreasing out of R, for 
instance, if the decreasing is exponential p oc exp{—r/rd) we 
can use as I. Then the cusp at v,nax transforms into a 
smooth peak. Since the free fall acceleration at r = ii is 
g = GM/R^, we can easily estimate the width of the peak: 
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(14) 



Near the Solar System Sv ~ (l/R) ■ 50 km/s. It is natural to 
assume that I <^ R (actually, it follows from supposition 4 
in the Introduction), so the distribution is still very narrow. 

To illustrate the above reasoning, let us consider the 
following density profile: 



const. 



if r < _R; 



dM/dr^^ f r-R\ (15) 
1 const ■ exp | | , it r > it 



rd 

In this case, we cannot determine the distribution ^(r) ana- 
lytically. However, if rd <C R, we can neglect the change in 
the gravitational field on the scale of rd and then find an ap- 
proximate solution. Fig. [T] represents distributions (|13|) (for 
the sharp halo edge, solid line) and for density profile (|15p 
(dotted line, we accepted rd = O.IR — 21 kpc). One can see 
that in the case of smoothed halo edge the peculiarity at 
Vmax — 562 km/s is transformed into a smooth peak, which 
is, however, quite narrow, and the general shape of the dis- 
tribution remains the same. We will discuss the influence of 
the halo density profile on the velocity distribution in the 
beginning of the Discussion section. 

The angular momentum of the dark matter particles 
hardly can be exactly equal to zero. If a particle possesses 
some specific momentum p, its velocity in gravitational field 
((2)1 is equal to: 



Vp = -; 
r 



Vr = 



(16) 



Since H16p is distinct from ((S)), distribution (jlip . strictly 
speaking, is not valid anymore. However, we may use it if 
the difference between (|16[) and ((5]) is small. Mathematically 
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Figure 2. The ratio between the double amphtudes 2A = Imax — 
Imin of direct detection signals II19I I calculated for anisotropic 
distribution ijlSp and Maxwell distribution p). 



it can be written as: 
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As we have already discussed in the Introduction, specific 
angular momentum of the majority of the dark matter par- 
ticles hardly can be larger than 900 kpc • km/s. Substitut- 
ing ro = R, r = tq to the inequality, we can see that 
its right part is equal to 110 km /s , while its left part 
is ~ 550^ km^/s^. Hence the inequality asserts near the 
Solar System, the influence of the angular momentum on 
the radial dynamics is still negligible for the majority of 
DMPs, and we can use (|lip up to rQ as before. Thus the 
DMP distribution throughout Vr in this approximation co- 
incides with (|13|) . However, the particles have also some dis- 
tribution throughout Vp. For simplicity we will suppose that 
/ oc exp{—Vp/2aQ) where ao = 80 km/s, though the distri- 
bution can be much narrower. Then the normalized DMP 
distribution near the Solar System can be closely approxi- 
mated by 



exp 



f{v) = 



2ag 



(18) 



where Vr G [—Vmax; Vmax], Vmax — 562 km/s. Distribution 
(|18p is strongly anisotropic and actually describes two col- 
liding beams of particles. 



3 DISCUSSION 

Fig.[T]represents distributions (|13|) and ([l]) (solid and dashed 
lines, respectively). One can see that (|13|) is much narrower 
and has much higher average velocity. The physical reason 
of it is obvious: in the case of Maxwell distribution Q the 
particles move almost circularly, which is why only a few of 
DMPs from the edge of the halo reach the Solar orbit. On 
the contrary, in the case considered in this article the ma- 
jority of DMPs comes from the halo edge and thus are much 



more accelerated by the gravitational field. Consequently, 
the question of what of the distributions, IT} or (|13|) . is cor- 
rect, can be reduced to whether the particles from the halo 
edge can reach the Solar orbit or not. In addition to the 
arguments presented in the Introduction we note that, ac- 
cording to H16p . a particle falling from r — R should have 
a specific angular momentum jj, ~ 4000 kpc ■ km/s, lest the 
particle can reach r — 8 kpc. This value is huge, it far ex- 
ceed not only the characteristic momentum of halo objects, 
but even the momentum of the disk, and thus looks very 
unlikely. So particles from the edge of the halo freely reach 
the Earth, and their spectrum should be closer to p3p . 

A similar consideration allows us to examine the depen- 
dence of velocity distribution (|13p on the density profile. Our 
assumption of the existence of a large region with p oc r"^ 
is approximately correct for massive spiral galaxies, such as 
the Milky- Way (Dutton et al,; 201Q). However, we obtained 
113} on the additional assumption that the edge of the halo 
is more or less sharp. Meanwhile, the outer region of the halo 
can have a density profile steeper than r~^, but not steep 
enough to be considered as a cutof f. As an instance, o ne can 
consider a double power- law halo (|Lisanti et al.|[201ll ). How 
can it influence on the velocity profile? The answer depends 
on the mass fraction of this steeper region with respect to 
the total mass of the halo. If the fraction is small, the dis- 
tribution differs little from (|13|l . as we demonstrated with 
distribution p5|l . Let us consider the case when the fraction 
is significant. We indicate the radius where the profile gets 
steeper than r~^ by 9\; hereafter we will name 'outer halo' 
the region out of JH. As we could see, in the model with 
sharp halo edge the majority of the particles comes from 
the edge of the halo. Expressing this fact mathematically, 
distribution ^(ro) is small for ro < R and goes to infinity at 
the edge of the halo (|ll|l . It is easy to show that in the case 
of the presence of a massive outer halo dark matter parti- 
cles mainly come from it, and the fraction of the particles 
with ro < is respectively small. Let us consider a system 
of particles with ro > 5H. According to ((Tjl, their contribu- 
tion to the halo mass in interval dr depends on r only as 
v~^{r). v~^{r), however, changes rather slowly inside the 
region where p oc r~^, since the potential there (see ([2])) 
depends on r only logarithmically. Therefore, the particles 
falling from the edge of halo provide almost the same con- 
tribution to the halo mass on each radius inside the region 
dM ~ const, which corresponds to p ~ r~^. Function ^(ro) 
should be chosen so that it reproduces the density profile, in 
particular, it should provide p oc r~^. However, as we could 
see, the particles from the outer halo by themselves give a 
very similar profile, and we need relatively few of particles 
with ro < in order to make it exactly r~^ . Consequently, 
^(ro) is small for ro < $H, and a significant fraction still 
comes from the halo edge. Thus, this property of the distri- 
bution does not depend on the exact density profile, being 
only a result of the assumption of strong anisotropy of the 
velocity distribution and of the flatness of potential ((2]). 

Function ^(ro) unambiguously determines the velocity 
distribution of the dark matter particles, and the above- 
mentioned common properties of ^(ro) directly correspond 
to characteristics of f (v). For the M i lky W ay galaxy is 
large {9\ > 100 kpc, iDutton et all (|2010l ')). which corre- 
sponds to w ~ 500 km/s for a terrestrial observer. Below 
this speed the distribution as a whole is similar to ((13} , and 
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it does not depend much on the density profile of the outer 
halo, since it is created by the particles with ro < fH. On 
the contrary, the distribution in region Vesc > v > 500 km/s 
strongly depends on the density distribution on the edge of 
the halo and can differ drastically from (|13|l . However, now 
we can estimate what distribution f(v) looks hke if profile 
p cc is valid only up to ~ 100 kpc, and then the 
halo has a massive outer region, for instance, a Navarro- 
Frank- White tail. We can expect that f{v) is similar to (|13p 
below V — 500 km/s and is totaly defined by the density 
profile of the outer halo for Vesc > v > 500 km/s. Since the 
outer halo is rather extensive, we can expect that a cusp in 
()13p is strongly smoothed. Fig. [1] illustrates all these prop- 
erties, though we made rather a small modification of the 
density profile H15[) . So the main characteristic features of 
distribution (|13|l are not sensitive to the density profile: the 
distribution is completely not Maxwell, a significant frac- 
tion of particles comes from the edge of the halo forming a 
high-velocity bump in f{v). 

The difference between U]) and (|13|) is important for 
various aspects of the dark matter physics. In the case 
of direct dark matter search the signal, roughly speaking, 
can be represented as a product of a part depending al- 
most not at all on the DMP dis tribution and an integral 
(|Belanger. Nezri. fc Pukhovll2009l ) 

/(„)= n iMdv (19) 

Here Vmin is the minimal DMP speed, to which the detector 
is sensitive, f{v) is the distribution in the Earth's frame of 
reference, obtained fr om ifH or l|13[l by a Gali lean transfor- 
mation (see details in lCerdefio fc GreenI (|2010l ). section 3.3). 
Because of the Earth's orbital motion I varies with a year 
period, and it is this variation that is observed in the direct 
detection experiments. Fig. [2] shows the ratio between the 
double amplitudes 2A = Imax — Imin of direct detection sig- 
nals calculated for anisotropic distribution (jlSp and Maxwell 
distribution ([1]), as a function of v,nin- One can see a very 
significant difference. 

The difference between distributions ^ and (I13p can 
also be important for the indirect dark matter search. Neu- 
trino observations, for instance, are trying to detect the dark 
matter annihilation going on in the centre of the Sun. The 
signal depends on the number of the DMPs captured by the 
Sun. Very roughly speaking, it is proportional to /(O). One 
can see that for distribution lfT3|l /(O) is approximately 20 
times smaller than for the conventional Maxwell distribution 

If the s-channel of the dark matter annihilation domi- 
nates (which is typical), (av) ~ const, and the signal is not 
sensitive to DMP distribution. However, if the p-channel pre- 
vails, a ~ const, and the signal is proportional to the aver- 
aged velocity of the particle collision, (vc) — O.Svmax in the 
case of distribution (|13|l and (vc) — x/S/tt Vorb in the case of 
Maxwell distribution HJ. Vmax relates to Vorb by dTOjl . and 
one can see that (vc) predicted by anisotropic distribution 
(|13p is much lower at the halo edge and much higher in the 
central region than the Maxwell one. Near the Solar System, 
however, they are nearly equal. Finally, Sommerfeld effect is 
inversely proportional to (vc) , and if it plays any role in the 
dark matter annihilation, it is also sensitive to the particle 
distribution. 
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